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Abstract 

A complete one loop analysis of the chromoelectric dipole contribution to the 
electric dipole moment (edm) of the quarks and of the neutron in N=l super- 
gravity including the gluino, chargino and neutralino exchange contributions and 
exhibiting the dependence on the two CP violating phases allowed in the soft 
SUSY breaking sector of minimal supergravity models is given. It is found that 
in significant parts of the supergravity parameter space the chromoelectric dipole 
contribution to the neutron edm is comparable to the contribution from the electric 
dipole term and can exceed the contribution of the electric dipole term in certain 
regions of the parameter space. An analysis of the contribution of Weinberg's 
purely gluonic CP violating dimension six operator within supergravity unifica- 
tion is also given. It is found that this contribution can also be comparable to the 
contribution from the electric dipole term in certain regions of the supergravity 
parameter space. 



1 Introduction 



The electric dipole moment (edm) of fermions is one of the important windows 
to new physics beyond the Standard Model (SM). In the SM the edm for the 
fundamental fermions arising from the Kobayashi-Maskawa CP violating phase 
is much smaller[0] than the current experimental limit of 1.1 x 10~ 25 ecm[]2| and 
beyond the reach of experiment in the forseeable future. In supersymmetric unified 
models new sources of CP violation arise from the complex phases of the soft 
SUSY breaking parameters which contribute to the edm of the quarks and of the 
leptons[3-9]. For the case of the quarks and the neutron there are also the color 



dipole operator and the CP violating purely gluonic dimension six operator [10, 11 



which contribute. However, with the exception of the work of ref. [[| virtually all 
previous analyses of the neutron edm have been done neglecting the contributions 
of these additional operators. The reason for this neglect is the presumption|| 
|7], || that their contribution to the neutron edm is small. In the analysis of this 
paper we show that the contributions of the color dipole operator and of the 
purely gluonic operator are not necessarily small and over a significant part of the 
supergravity parameter space their contribution to the neutron edm is comparable 
to the contribution from the electric dipole operator and can even exceed it in 
certain regions of the parameter space. 

In this Letter we first derive the full one loop contribution to the color dipole 
operator arising from the gluino, the chargino and the neutralino exchange con- 
tributions exhibiting the dependence on the two CP violating phases allowed in 
the soft SUSY breaking sector of minimal supergravity models. To our knowledge 
this is the first complete one loop analysis of this operator. We also recompute the 
purely gluonic dimension six operator to incorporate the dependence on the two 
CP violating phases. We then give a numerical analysis of the relative strength 
of the color dipole contribution to the neutron edm relative to the contribution 
from the electric dipole term using the framework of supergravity unification under 
the constraint of radiative breaking of the electro-weak symmetry. Our analysis 
shows that the color dipole contribution relative to the electric dipole contribu- 
tion can vary greatly over the parameter space of the model. We find that there 
are significant regions of the parameter space where the color dipole contribution 
can be comparable to the contribution from the electric dipole term and in some 
regions of the parameter space the color dipole contribution can even exceed the 
contribution from the electric dipole term. A similar analysis is also carried out 
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for the purely gluonic dimension six operator. Here again one finds that there 
exist regions of the parameter space where the contribution of the purely gluonic 
dimension six operator to the neutron edm can be comparable to, and may even 
exceed, the contribution from the electric dipole term. Thus one is not justified 
in discarding the effects of the color dipole operator and of the gluonic dimension 
six operator in the neutron edm computation. 

The outline of the Letter is as follows: In Sec. 2 we give the complete one 
loop calculation including the gluino, the chargino and the neutralino exchange 
contribution to the color dipole operator and discuss its contribution to the electric 
dipole moment of the neutron including its dependence on the two CP violating 
phases in minimal supergravity. In Sec. 3 we recompute the dimension six gluonic 
operator to include the effects of the two CP violating phases. In Sec. 4 we give the 
renormalization group analysis in minimal supergravity of the relative numerical 
strengths of these contributions and show that they can be comparable to the 
electric dipole contribution in significant regions of the parameter space. 

2 Analysis of Chromoelectric Dipole Contribu- 
tion 

The parameters that enter in minimal supergravity grand unification with radiative 
breaking of the electro- weak symmetry can be taken to be mo, mi/2, Aq, tanft and 
phase (/x), where m is the universal scalar mass, mw 2 is the universal gaugino mass, 
A is the universal trilinear coupling, tan/? =< H 2 > / < Hi >, where H 2 gives 
mass to the up quark and H\ gives mass to the down quark, and /i is the Higgs 
mixing parameter |12], [l~3| . As noted in Sec. 1 only two phases in the soft SUSY 
breaking sector have a physical meaning in minimal supergravity, and we choose 
them to be the phase aAo of A and the phase # M o of /Jo- The quark chromoelectric 
dipole moment is defined to be the factor d c in the effective operator: 

Cj = - l -d c qa, ul ,T a qG^ a (1) 

where T a are the generators of SU(3). In the following we give an analysis of 
the one loop contributions to d c from the chargino, the neutralino and the gluino 
exchange diagrams shown in Figs, la and lb. 
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2.1 Gluino Contribution 

The quark-squark-gluino vertex can be derived using the interaction ]T5 



- C q ^- g = V2g s T« k {-ql l -^M- R + q^^M- L ) + H.c, (2) 

i=u,d 

where a = 1 — 8 are the gluino color indices, and j, k — 1 — 3 are the quark and 
squark color indices. The complex phases enter via the squark (mass) 2 matrix M| 
which can be diagonalized by the transformation 

D\M 2 D q = ^g(M 2 1) M 2 2 ) (3) 

where 

<?l = Dqiiqi + D q i 2 q 2 (4) 

Qr = D q21 qi + D q22 q 2 . (5) 

and where q\ and q 2 are the mass eigenstates. Writing C in terms of q~\ and q 2 and 
integrating out the gluino and squark fields and by using the identities 

and 



^ -[S u 6 kj - -6M (6) 



3 

jabcrpbrpc _^~£ a j^j) 



one can obtain the gluino exchange contribution to d c . We find 



2 ^2 

d q -gluino = -J-T Yl Im ( r f)jjl C ( mI")' ^ 

k=l <lk Qk 

where r* fc = D q2k D* lk and rrig is the gluino mass and 

. . 1 , 2rlnr 18/nr\ . , 

6(r — 1) 1 — r 1 — r 

2.2 Neutralino Contribution 

Here the CP violating phases enter via the the squark (mass) 2 matrix, which 
contains the phases of A q and n, and via the neutralino mass matrix given by 



/rhi — M z sin 9w cos/5 M z sin 9w sin (3 \ 

m 2 M z cos 6w cos /3 — M z cos 6^ sin (3 

— M z sm9w cos/5 M z cos 6*^ cos /3 — fi 

V M z sin 9w sin /3 — M z cos^iy sin/3 — // / 

(10) 
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which carries the phase of fi. The matrix M x o is a complex non hermitian and 
symmetric matrix, which can be diagonalized by a unitary transformation such 
that 

X T M x oX = diag(m x o, m x o, m x o, m x o) (11) 



Quark-squark-neutralino vertex can be derived from the interaction |13 

4 
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+((3 uj D ull - 5 uj D u21 ) 1 + ^ 5 ]x°jUi+V2u[(a uj D ul2 - A.22) 1 ^ 

+ (/3ujD ul2 - 5 U3 D u22 )^^]x]u 2 + (u^d) + H.c. (12) 
where a, (3, 7 and 5 are given by 

2mv^ sin j3 (cos /5) 



= „;„ /3/„„„ m l 1 ^ 



= eQ„( d )X 1 * H — X 2 *(T 3u ( d ) — Q u (d) sin 2 6^) (14) 

COS ts'W 



gQ u{ d) sin 2 6 W 
cos 9 w 

K W = . 9m ^-' (16) 

2m w sin /3 (cos /3) 



lu{d)j = eQ u{d) X 1:j — X 2j (15) 

tub 



and where 



X' Xj = X X j cos 6 W + X 2j sin 6 W (17) 
X' 2j = -Xij sin 8 W + X 2 j cos 8 W . (18) 
The one loop analysis using Eq.(12) gives for d c 

dq-neutralino = £ £ M^fe) —J) , (19) 

iD7r fe=li=l 1V1 qk 



where 



Here 



rj gik = [-V2{tan 6 w (Q q - T 3q )X u + T 3q X 2i }D* lk + K q X bi D* 2k \ 

(V2tan9 w Q q XiiD q2 k - K q X hi D qlk ). (20) 



m u rrid 
\'2m w sinf] \ / 2m w cos/3 



where 6 = 3(4) for T 3q = — ~(|) and : 



1 Irlnr . , , 
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2.3 Chargino Contribution 

Here the CP violating phases enter via the squark (mass) 2 matrix and via the 
chargino mass matrix given by 

(_ m 2 V2m w smp\ 

W cos p /i / 

which involves the phase of //. The chargino matrix can be diagonalized by the 
unitary transformation: 

U*M c V~ l = diag(m x +, m x +) (24) 



The quark-squark-chargino vertex can be derived from the interaction [13 

1+75 



- Cq-q-x+ = gU^UuDdu - K d U 12 D d21 ) 



2 



-(K u V* 2 D du )^—^-]xfdi+gu[(U 2 iD du - K d U 22 D d21 ) - ^ 7s 
-K u V 22 D dll )^—^-}xidi+gu[(U u D dl2 - n d U 12 D d22 ) 1 + 75 
-(tt u V* 2 D dl2 ) 1 ^ }xtd2+gu[{U 2 iD dl2 - K d U 22 D d22 ) + ryb 
-(n u V 22 D dl2 ) 2 75 ] J 2 

+(« " d, ^ " V, X, + -> X, C ) + #c, (25) 
and for d c one gets for the up and down flavors 

du-chargino = "Tp-^ 12 12 M 1 **) J^T B ( wf") > ( 26 ) 
fc=l i=l Jfc dfc 

-g 2 g s 2 2 m x + m 2 + 

d d-char gl no = "TcZl" 12 12 Im ( F ^fc) T7+" B ( Jjt) > ( 27 ) 

iD7r fe=li=l iW «fe 

where 

r uifc = WuV^-DdifcCC^i-DSifc - KdU* 2 D d2k ) (28) 



r„ fc = ^Au^i^i* - K u V*D* u2k ). (29) 
The contribution to EDM of the quarks can be computed using the naive dimen- 
sional analysis |14[ which gives 

d « = (3°) 
where if is the renormalization group evolution of the operator of Eq.(l) from the 
electroweak scale down to hadronic scale §]. For the neutron electric dipole 
moment d n we use the naive quark model d n = (4d d — d u )/3. 
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3 CP Violating Purely Gluonic Dimension 6 op- 
erator 



The gluonic dipole moment d G is denned to be the factor in the effective operator 



£■1 — ~—d G fa/3^G a ^ p Gp u G 1 x a e fil ' Xa (31) 



where G afll/ is the gluon field strength tensor, f a/3l are the Gell-Mann coefficients, 
and e^ uXu is the totally antisymmetric tensor with e 0123 = +1. 

Dai et al [[□]] have calculated d G considering the top quark-squark loop with 
a gluino exchange in terms of the complex phase <fi of the gluino mass while the 
squark (mass) 2 matrix in their analysis was considered to be real. In our case the 
squark (mass) 2 matrix is complex and carries the phases of A q and of fi. We have 
recalculated d for our case and we find 

d G = -3a s m t (^) 3 Im(r t 12 )^^H(^,^ 2 ,z t ) (32) 



where 



and 



'47T m~ 



,M ? 2 _ f m t 



z t = O 2 (33) 



m- g m g 



V] 2 = D t22 D; i2 (34) 

where D t is the diagonalizing matrix for the stop (mass) 2 matrix, and the function 
H is the same as in |Tl| . The contribution to d n from d G can be estimated by the 
naive dimensional analysis which gives 

d£ = (35) 

where M is the chiral symmetry breaking scale with the numerical value 1.19 GeV, 
and r] G is the renormalization group evolution of the operator of Eq.(31) from the 
electroweak scale down to the hadronic scale, rf and rf" have been estimated to 



be ~ 3.3 15 



4 RG Analysis and Results 

We discuss now the relative size of the three contributions to the neutron edm, 
i.e., from the electric dipole operator, from the color dipole operator and from the 
purely gluonic operator. To study their relative contributions we consider a given 



6 



point in the supergravity parameter space characterized by the set of six quantities 
at the GUT scale: mo, mi , Aq, tan/?, 6^0 and a A o- 111 the numerical analysis we 
evolve the gauge coupling constants, the Yukawa couplings, magnitudes of the 
soft SUSY breaking parameters, fi and the CP violating phases from the GUT 
scale down to the Z boson scale. We use one-loop renormalization group equations 
(RGEs) for the soft SUSY breaking parameters and two-loop RGEs for the Yukawa 
and gauge couplings. Using the data gotten from the RG analysis at the scale Mz 
we compute the contributions to the quark edm from the electric dipole part {d^), 
from the color dipole part (d^), and from the purely gluonic part (d^)- These are 
further evolved to the hadronic scale by using renormalization group analysis as 
discussed in Sees. 2 and 3. 



Table 1: m~ g = 500 GeV, m = 2000 GeV, \A \ = 1.0 



case 


tan /3 


phases (rad) 


<(10" 26 ecm) 


c£'(10- 26 ecm) 


d% (10- 2e ecm) 


(i) 


2 


MO = 0.2, a Ao = -0.5 


0.124 


3.049 


-1.168 


(ii) 


2 


9^ = 0.2, a Ao = 0.5 


-8.46 


12.84 


17.315 


(iii) 


4 


M = 0.2, a Ao = -0.5 


-4.33 


0.764 


-5.13 


(iv) 


4 


6^ = 0.2, a Ao = 0.5 


-11.74 


-12.65 


7.28 



To exhibit the importance of the color dipole operator and of the purely glu- 
onic operator we display the relative sizes of d„, d% and of d^ for few illustrative 
examples in Table 1. One can understand the smallness of <if for case(i) in the 
following way. The main contributions to d% arise from the chargino and from 
the gluino exchange. The chargino exchange gives a negative contribution while 
the gluino exchange gives a positive contribution and the smallness of d„ is due 
to a cancellation between these two. For the color dipole term d% there is also a 
cancellation but this time the cancellation is only partial and it occurs between 
the d-quark and the u-quark contributions. Because of a large cancellation for d% 
and only a partial cancellation for d% one has dominance of \d%\ over \d%\ in this 
case. 

One may contrast the result of case(i) with that of case(ii) where the sign of 
a Ao is switched. Here the gluino contribution in d% switches sign and this time one 
has a reinforcement of the chargino and the gluino contributions making \d^ \ much 
larger than for case(i). Further, in the color dipole part the d-quark contribution 
in the gluino exchange switches sign and there is a reinforcement of the d-quark 
and the u-quark contributions making \d%\ larger than for case(i). We see then 
that in this case \d%\ and \d%\ are comparable. One may note that a very large 
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change occurs for the purely gluonic term in going from case(i) to case(ii). To 
understand the large shift in the value of we display explicitly the imaginary 
part of Eq. (34) 

Im ( r t 12 ) = / M? ^* 2 r (m \A t \ sin a A + sin 6^ cot (3), (36) 
v ti iv 

where 9^ and a a are the values of 9^ and of aA at the electro-weak scale. From 
Eq. (36) we see that the magnitude of Im(r t 12 ) depends on the relative sign and 
the magnitude of 9^ and of a a- Thus one has a cancellation between the A t and 
the \i terms when 9^ and a a have opposite signs as in case(i) and a reinforce- 
ment between the A t term and the \i term when 9^ and aA have the same sign as 
in case(ii). Thus one can qualitatively understand the largeness of d^ in case(ii) 
relative to case(i). The very large reduction in d% for case(iii) relative to case 

(ii) occurs because of a reduction in the down quark contribution, which is the 
dominant term in d n , due to a switch in the o:a sign and a change in the value of 
tan/3. In going from case(iii) to case(iv) the further switch in the sign of aA Q once 
again increases the down quark contribution making d% the largest in magnitude. 
It is interesting to note that each of the four cases of Table 1 leads to a distinct 
pattern of hierarchy among the three contributions. Thus one has 

(i) ICI > 1^1 > Kl (h) \d°\ > \d c n \ > K\ 

(iii) 1^1 > K\ > \d% (iv)K| > \d*\ > |d*| 

The analysis clearly shows that any one of the three contributions can dominate 
d n depending on the part of the parameter space one is in. We note also that d% 
for all the four cases in Table 1 is consistent with experiment while the total edm 
which includes the color and the purely gluonic part for cases (ii) and (iv) would 
be outside the experimental bound. 

In Fig. 2 we display the magnitudes of df t , d% and d% as a function of rrig 
for a specific set of m , A , tan/3, aA and 9^ values. Here we find that \d%\ is 
comparable to \d%\ over most of the rrig region and in fact exceeds it for values of 
rrig below ~ 800 GeV. Further \d%\ also exceeds \d%\ for values of m g below ~ 400 
GeV in this case. The broad peak in d^ arises from a destructive interference 
between the gluino and the chargino components of d^ which leads to a relatively 
rapid fall off of <if for values of m g below ~ 500 GeV. There are other regions where 
similar phenomena occur. To exhibit the commonality of the largeness of d% we 
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give in Fig. (3a) a scatter plot of the ratio \d^/d^\ for the ranges indicated in the 
Fig. (3a) caption. We see that there exist significant regions of the parameter space 
where the ratio \d^/d^\ ~ 0(1) and hence d% is non-negligible in these regions. A 
scatter plot of \d^/d^\ is given in Fig. (3b) and it exhibits a similar phenomenon. 

In conclusion we have given the first complete one loop analysis of the chro- 
moelectric contribution to the electric dipole moment of the quarks and of the 
neutron exhibiting the dependence on the two arbitrary CP violating phases al- 
lowed in minimal supergravity unification. We find that the relative strength of 
the chromoelectric contribution varies sharply depending on what part of the su- 
pergravity parameter space one is in. In significant parts of the parameter space 
the chromoelectric contribution is comparable to the electric dipole contribution 
and can even exceed it in certain regions of the parameter space. A similar con- 
clusion also holds for the contribution of Weinberg's CP violating purely gluonic 
dimension six operator. The analysis of the neutron edm exploring the full param- 
eter space of supergravity unified models is outside the scope of this Letter and 
will be discussed elsewhere. We only state here that the inclusion of all the three 
components of the neutron edm, i.e, d^ , d%, and d%, is essential in making reliable 
predictions of the neutron edm. In fact, the full analysis shows there exist regions 
of the supergravity parameter space where internal cancellation among the three 
components lead to an acceptable value of the neutron edm without either the 
use of an excessively heavy SUSY spectrum or an excessive finetuning of phases. 
These results have important implications for the mechanisms needed to suppress 
the neutron edm in SUSY, for the effect of CP violating phases on dark matter 
and on the analyses of baryon asymmetry in the universe. 
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6 Figure Captions 

Fig. la: One loop diagram contributing to the color dipole operator where the 
external gluon line ends on an exchanged squark line in the loop.Squarks are rep- 
resented by qk in the internal lines. 
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Fig. lb: One loop diagram contributing to the color dipole operator where the 
external gluon line ends on an exchanged gluino line labelled by g in the loop. 
Fig. 2: Plot of d„, d%, and as a function of m g when m =2000 GeV, tan/3=3.0, 
|A)| = l-0, # W) =0.1 and aA =0-5. 

Fig. 3a: Scatter plot of the ratio \d^/d^\ as a function of mo for the case |A | = 1.0, 
tan/3 = 2.0 and the other parameters in the range 200 GeV< m g < 600 GeV and 
-7r/5 < 9 w ,a Ao < 7r/5. 

Fig. 3b: Scatter plot of the ratio \d^/d^\ as a function of mo for the same range 
of parameters as in Fig. 3a. 
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